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ABSTRACT 

We present VLT/FORS2 time-series spectroscopy of the Wolf-Ray et star #41 in the Sculptor 
group galaxy NGC 300. We confirm a physical association with NGC 300 X-1, since radial 
velocity variations of the He II A4686 line indicate an orbital period of 32.3 ± 0.2 hr which 
agrees at the 2a level with the X-ray period from Carpano et al. We measure a radial velocity 
semi-amplitude of 267±8 kms^^, from which a mass function of 2.6 ± 0.3 Mq is obtained. A 
revised spectroscopic mass for the WN-type companion of 26^5 Mq yields a black hole mass 
of 20 ± 4 A/0 for a preferred inclination of 60 — 75°. If the WR star provides half of the mea- 
sured visual continuum flux, a reduced WR (black hole) mass of 15^2 5 -^0 (14.5l2 5 -^o) 
would be inferred. As such, #41/NGC 300 X-1 represents only the second extragalactic Wolf- 
Rayet plus black-hole binary system, after IC 10 X-1. In addition, the compact object respon- 
sible for NGC 300 X-1 is the second highest stellar-mass black hole known to date, exceeded 
only by IClOX-1. 

Key words: (galaxies:) individual: NGC 300 - Stars: Wolf-Rayet - X-rays: binaries - X-rays: 
individual: NGC 300 X-1 



; 1 INTRODUCTION 

] High mass X-ray binaries (HMXB) typically comprise OB stars 
I plus either a neutron star or black hole, in which high X-ray lu- 
' minosities (~10'*** ergs^^) arise from accretion disks around the 
, compact object. Accretion disks are fed through a combination of 
■ Roche lobe overflow and stellar winds from the early-type compan- 
ion. Very few known HMXB are known to host black-holes, Cyg 
X-1 in the Milky Way, X-1 and X-3 in the Large Magellanic Cloud 
and X~7 in the Local Group galaxy M33. If a merger is avoided, 
the OB companion will potentially evolve through to a Wolf-Rayet 
phase, producing a system comprising a helium star plus a black 
hole or neutron star (Tutukov & Yungelson 1973). Indeed, van den 
Heuvel & de Loore (1973) proposed the Galactic HMXB Cyg X-3 
as a helium-star plus compact object system, which was observa- 
tionally confirmed by van Kerkwijk et al. (1992). Such systems are 
believed to be very rare, with as few as ~ 100 helium star plus black 
hole pairs in the Galaxy (Ergma & Yungelson 1998). 

To date, the only confirmed Wolf-Rayet plus black hole bi- 
nary system is IC 10 X-1 in the dwarf irregular galaxy IC 10. This 
system has a period of 34.9 hr, hosts a WN-type Wolf-Rayet star 
[MAC 92] 17A (Crowther et al. 2003) and an unseen companion 
which is currently the record holder amongst stellar mass black 



* Based on observations made with ESO telescopes at the Paranal Obser- 
vatory under programme ID 384.D-0093(A) 
f Paul.Crowther@shef.ac.uk 



Table 1. Log of VLT/FORS2 spectroscopic observations of #41 in 
NGC 300. UT dates and MJD's refer to the start of the 1,535 s exposures. 
We include individual radial velocities, Vr , as measured from Gaussian fits 
to He II A4686. Phases adopt a period of 32.3 hr, where phase refers to 
MJD55118.97559±0.01554). 



UT Date 


MJD 


DIMM 


sec z 


A4686 


Phase 




-55119 


ff 




kms^^ 




01:06 15 Oct 09 


0.046180 


1.2 


1.29 


347.0 ± 17.7 


0.052 


05:28 15 Oct 09 


0.227922 


2.3 


1.08 


364.6 ±25.5 


0.187 


01:40 16 Oct 09 


1.069495 


1.1 


1.18 


-64.2 ±16.9 


0.812 


05:58 16 Oct 09 


1.249073 


1.7 


1.14 


136.0 ±25.0 


0.945 


04:07 19 Oct 09 


4.171590 


0.7 


1.03 


369.0 ± 14.3 


0.114 


08:24 19 Oct 09 


4.350382 


0.8 


1.98 


456.6 ±20.1 


0.247 


02:31 20 Oct 09 


5.104896 


0.9 


1.06 


-68.3 ± 12.1 


0.807 


06:41 20 Oct 09 


5.278856 


0.7 


1.31 


109.4 ± 17.3 


0.936 


00:14 25 Oct 09 


10.010318 


1.2 


1.34 


325.3 ±25.7 


0.447 


06:28 25 Oct 09 


10.269520 


0.7 


1.33 


10.4 ± 14.3 


0.639 



holes, exceeding 23.1 ±2.1 Mq (Prestwich et al. 2007; Silverman 
& Filippenko 2008). 

A second candidate extragalactic black hole plus Wolf-Rayet 
system is NGC 300 X-1 (Carpano et al. 2007a), in the southern 
Sculptor group spiral galaxy NGC 300 which lies at a distance of 
1.88 Mpc (Gieren et al. 2005). This source is spatially coincident 
with #41 from the Wolf-Rayet catalogue of Schild et al. (2003) 
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Figure 1. Representative VLT/FORS2 spectroscopy of #41 revealing radial 
velocity variations in He II A4686 (gaussian fits are shown as dotted lines), 
successively offset by 1.5x10"^'^ erg s~^ cm~^ for clarity. 



which was confirmed as a WN-type Wolf-Rayet star by Crowther et 
al. (2007). Both systems exhibit similar X-ray properties (Carpano 
et al. 2007b; Barnard, Clark & Kolb 2008). However, no evidence 
exists to date of a physical link between NGC 300 X-1 and the 
Wolf-Rayet star #41. This is the purpose of the present Letter. 

New Very Large Telescope (VLT) optical spectroscopic time- 
series of #41 are discussed in §|2]in which variations are revealed 
in the Hell A4686 emission line. §[3]compares the inferred orbital 
period with X-ray lightcurves and derives a semi-amplitude for the 
Wolf-Rayet star, from which a mass function is obtained. §|4] pro- 
vides a revised mass for the Wolf-Rayet star, placing strict limits 
upon the mass of the compact companion. We conclude with a brief 
discussion in §|5] 



2 OBSERVATIONS 

Here we present new VLT optical spectroscopy of #41 obtained 
with the Focal Reducer/Low Dispersion Spectrograph #2 (FORS2) 
in multi-object spectroscopy (MOS) mode from 14—24 October 
2009. Two 1,535 s exposures were obtained on each of five non- 
consecutive nights using the 600B grism, centred at 465 nm. A log 
of our observations is presented in Table [T] including DIMM seeing 
measurements, which indicated typical seeing of 0.7-1.2 arcsec, 
except at two epochs for which the seeing exceeded 1.5 arcsec. 

MOS allows 19 sources to be simultaneously observed. We 
included 7 Hll regions in NGC 300 plus 12 WR candidates from 
Schild et al. (2003). Of these, 3 sources have previously been spec- 
troscopically confirmed as WR stars, namely #9, #41 and source 12 
from Bresolin et al. (2009). Other sources will be discussed else- 
where. 

After bias subtraction and flat field correction, a standard ex- 
traction was performed using IRAF, with wavelength and flux cali- 
bration carried out using FIGARO. 1.0 arcsec slits provided a spec- 




0.25 0,50 0,75 1,00 1,25 1,50 

J^■requency (cycles/day) 

Figure 2. String-length periodogram of the VLT/FORS2 radial velocities. 
The minimum string length indicates the best period, at a frequency of 
0.7422 cycles/day. The solid, dashed and dotted lines represent the 68.3%, 
95.4% and 99.7% confidence levels, respectively (see text for details). 



tral resolution of 4.2 A - as measured from comparison arc lines - 
with a wavelength coverage of 3300-5800A for #41. Faint H/3 and 
[Olll] AA4959, 5007 nebular lines were detected in the #41 slit. 
#41 does not appear to be responsible for these features since the 
intensity of the nebula is weak and spatially close to #41, and has 
a lower systemic velocity of ~95 kms^^ (versus 202 kms^^ for 
#41). 

Spectrophotometric standard stars were observed with F0RS2 
in long sht mode on the nights on 15 Oct (HD 49798), 18-19 Oct 
(LTT2415) and 24 Oct (LTT 7987), providing a wavelength cover- 
age of 3300-6230A. An absolute flux calibration was achieved for 
#41 using B = 22.71 mag from Carpano (2006). 



3 ORBITAL PERIOD 

Several VLT/F0RS2 spectroscopic observations in the vicinity of 
Hell A4686 are presented in Fig.[T] revealing large radial velocity 
variations. Gaussian profiles are fit to individual A4686 profiles, 
with individual centroids listed in Table[T] 

In view of the sparsely sampled datasets, we have employed 
the string-length approach of Dworetsky (1983). The data are 
folded on a set of trial frequencies and the total length of "string" 
required to join the observations in phase order is calculated. The 
smallest string length found from the search is assumed to corre- 
spond to the correct period. The resulting periodogram is shown in 
Fig-IH The deepest trough conesponds to a period of 32.3 ± 0.2 hr. 
The error on this period was computed by constructing 10 000 syn- 
thesized datasets and measuring the standard deviation of the po- 
sitions of the deepest troughs in the resulting periodograms. The 
synthesized datasets were obtained by "jiggling" each data point 
about its observed value by an amount given by its error bar multi- 
plied by a number output by a Gaussian random-number generator 
with zero mean and unit variance. 

It can be seen that our derived period is in agreement at the 
2a level with the Swift X-ray period of32.8±0.2hr(l(T; Carpano 
et al. 2007b), which gives us confidence that we have identified the 
correct value. Moreover, the minimum string length of our period 
(1.49) compares favourably with the string length of a perfect sinu- 
soid ( 1 .46) and the string length of a sinusoid with noise consistent 
with the error bars on the observed data added to it (1.48). 

To further test the significance of our derived period, we used a 
randomization technique (Fisher 1935). The radial velocities were 
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randomly reassigned to the times of observation, thiereby preserv- 
ing the data sampling and the mean and standard deviation of the 
original dataset. A set of 10000 randomised datasets were con- 
structed in this way and then subjected to the same string-length 
periodogram analysis. By constructing a cumulative distribution 
function of the resulting minimum string lengths, we are able to 
place confidence limits on the significance of a given string length, 
as shown by the horizontal lines in Fig. (2] We are able to reject the 
hypothesis that our preferred period is due to noise at 99.43% con- 
fidence and eliminate the next highest troughs as most likely due 
to noise. Note that the troughs around 0.25 and 1.75 cycles/day are 
due to the one-cycle-per-day alias. As a check on the string length 
method, we also computed a Lomb-Scargle periodogram (Press & 
Rybicki 1989) and obtained consistent results. 

Adopting our optical period, we present the phased radial ve- 
locity measurements in Fig. [3] with a semi-amplitude of K2 = 
267.5±7.7 kms"\ In the Fi gure, phase corresponds to MJD 
55118.97559±0.01554. 

Armed with the semi-amplitude and orbital period, we may 
now derive the system mass function, 

^, , PKi Misin^i 

where Mi is the compact object mass, and q = M2 /Mi . The de- 
rived mass function is /(m) — 2.63 ± 0.33 M©, and would corre- 
spond to the compact object mass in the case of a system viewed at 
an inclination of 90° whose companion mass is negligible. In the 
case of a massive companion star with q ^ 1, the compact object 
would have a minimum mass of 4 f{m). As such, the compact ob- 
ject in NGC 300 X-1 is indeed a black hole, such that this system 
represents only the second confirmed WR plus black hole binary. 



4 WOLF-RAYET PROPERTIES 

We present our new, combined (phase-corrected), rectified 
VLT/F0RS2 spectrum of #41 in Fig. |4] This high quality spec- 
trum confirms a weak-lined WN5 subtype, previously inferred by 
Crowther et al. (2007) from lower resolution, lower S/N spec- 
troscopy obtained with VLT/F0RS2 using the 300V grism in Jan- 
uary 2007. Overall, the visual spectrum of #41 is similar to other 
weak-lined WN5 stars, namely WR 49 in the Milky Way and Brey 
65b (= NGC 2044 West 5C) in the LMC, taken from Hamann et 
al. (1995) and Walborn et al. (1999), respectively. The He II A4686 
equivalent width (Wx ~56A), and line width (FWHM~17A) in 
#41 are somewhat lower than LMC and Milky Way counterparts, 
Wx=l 10-140A and FWHM=22-24A. 

In order to reassess the mass of #41 we have calculated a 
synthetic model using the Hillier & Miller (1998) line-blanketed, 
non-LTE model atmosphere code. With respect to Crowther et al. 
(2007), somewhat more sophisticated atomic models are consid- 
ered, namely H, He, C, N, O, Ne, Si, P, S, Ar, Fe and Ni. Elemen- 
tal abundances are set to 40% of the solar value (Urbaneja et al. 
2005), with the exception of H and CNO elements. Clumping is 
accounted for, albeit in an approximate manner, with a (maximum) 
volume filling factor of 10%, such that the derived mass-loss rate 
is three times smaller than the value that would have been obtained 
by assuming an homogeneous wind. 

In view of the weak He I line spectrum in #41, we have based 
our analysis upon Hell (A4686, A5411) and Niv-V (AA4603-20, 
A4058, A7103-7129) line diagnostics. Overall good agreement is 
found, which is remarkable in view of the close proximity of the 
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Figure 3. Radial velocity vaiiations of A4686 He II phased to 32.3 hr, from 
which a systemic velocity of Vr = 202 ± 7 kms^^ and semi-amplitude of 
K2 = 267.5 it 7.7 kms^^ is obtained. 
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Figure 4. Comparison between rectified, phase corrected VLT/FORS2 
spectrum of #41 with Galactic (WR 49) and LMC (Brey 65b = NGC 2044 
West 5C) weak-lined WN5 stars, respectively from Hamann et al. (1995) 
and Walborn et al. (1999). 



black hole to #41. The only significant discrepancies are that N III 
A4634-41 is not reproduced in the synthetic spectrum and excess 
emission is observed in the upper Pickering-B aimer series, the lat- 
ter potentially arising from the accretion disk. 

In Fig. |5] we present our new combined flux calibrated 
spectrum of #41, together with re-calibrated spectroscopy from 
Crowther et al. (2007) for A > 5800A. An optimum fit to the 
spectrum of #41 is included in the figure, and reveals the fol- 
lowing stellar parameters - T* ~ 65 kK, log(L/Lo) ~ 5.92, 
M ~ 5 X 10"*^ Mq yi~^, Vao ~ 1300 km s~^, plus a nitrogen mass 
fraction of ~0.5%, with negligible hydrogen adopted. With respect 
to Crowther et al. (2007), the main revision relates to a reduced 
absolute magnitude of My = -5.0 mag, on the basis of a lower 
interstellar reddening of E{B — V) = 0.4 mag. T, should be reli- 
able to ±5 kK, resulting in uncertainties of ±0.2 mag in bolometric 
corrections. Together with ±0.05 mag uncertainties in E{B — V), 
stellar luminosities should be reliable to ±0.14 dex. 

From our derived parameters, we obtain a spectroscopic WR 
mass of 26I5 ^0 basis of the Schaerer & Maeder (1992) 

mass-luminosity relation for hydrogen-free WR stars. The principal 
uncertainty in our inferred WR mass relates to the absolute visual 
magnitude of the WR star. Since our adopted visual magnitude is 
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Figure 5. Combined, phase-corrected VLT/FORS2 spectroscopy of #41 
from Oct 2009 and Jan 2007 (black) with synthetic spectrum overlaid (red), 
reddened by E(B — y)=0.4 mag, as described in §|4] 



based upon ground-based imaging, it is possible that other contin- 
uum sources are included in the photometry. In this case, the WR 
emission line equivalent width would be diluted by the continuum 
of other nearby sources. 

Indeed, while #41 closely resembles other weak-lined WN5 
stars (recall Fig.|4](, its He II A4686 equivalent width is indeed lower 
by a factor of ~2. This hints at a potential factor of two line dilu- 
tion from unresolved companions along the sight-line towards #41. 
In this case, the WR luminosity would be reduced to log(L/LQ) 
= 5.57±0. 14 with the mass-loss rate unaffected, implying a spec- 
troscopic mass of 15^2 fjMq. For reference, Hamann, Grafener, & 
Liermann (2006) estimated spectroscopic masses of 15-19 Mq for 
weak-lined WN5 stars in the Milky Way. In view of these issues, 
we shall evaluate black hole masses using values of both 15 and 26 
Mq for the WN star. 



5 DISCUSSION AND CONCLUSIONS 

In Table|2]we present black hole masses for inclinations of 45, 60 
and 90° for our favoured 26 Mq spectroscopic WR mass, plus the 
15 Mq case, resulting from the WR star contributing 50% of the 
visual continuum light. We include the separation between the com- 



Table 2. Derived black hole mass. Mi, in NGC 300 X-1 for i=45, 60° or 
90°, for cases in which the WN star contributes either 50% (M2 = 15Mq) 
or 100% (M2 = 26Mq) of the visual light. We include the separation 
between the components, a, and the WR radius, R2, as a fraction of the 
Roche lobe radius, rj^ (Eggleton 1983). 



M2 (WR) 

Mq 


R2 (WR) 

Rq 


i 


Ml (BH) 

Mq 


a 

Rq 


Ra/rL 


I5+* 
'-'-2.5 

15+" 
'-'-2.5 

'-'-2.5 


4.8±0.1 
4.8±0.1 
4.8±0.1 


45I 

60 

90^ 


21.5+1-^ 
15.6li-| 

n.6tli 


17.0+^7 
16 1+1° 
15 5+1.0 
'^■-•-0.8 


0.81 
0.80 
0.79 


26+1 
26+1 
26tl 


7.2±0.1 
7.2±0.1 
7.2±0.1 


451 

60 

90^ 


27.8t^j 
20.6+2? 
16.9t?i 


19.4i-g 


0.99 
0.98 
0.96 



1 The apparent glancing eclipse of the X-ray emitting accretion disk sug- 
gests j = 60 - 75°. 

ponents using Kepler's third law in each case, obtained from 



[G'(Mi + M2)]^'''' 
corresponding to only 2.5-3.5 Wolf-Rayet radii. The mass ac- 
cretion rate required to sustain Lx = 2 xlO''* ergs~^ is 3.5 
X 1O~*M0 yr~^, if the adopted efficiency of gravitational release 
is ~10% (see Shakura & Sunyaev 1973). This is <1% of our de- 
rived mass-loss rate of #41. However, Table |2] also shows that the 
WR star would completely fill its Roche lobe for the 26 Mq case 
(i <35° is excluded), and equates to 80% of its Roche lobe radius, 
r'L (Eggleton 1983), for the 15 Mq case. Therefore, the accretion 
disk may be fed primarily through Roche lobe overflow. For com- 
parison, the higher temperature obtained for the WN star in IC 10 
X-1 by Clark & Crowther (2004) would favour a wind-fed accre- 
tion disk, since the WR radius is ~ 0.5rL in that system. 

IC 10 X-1 is an eclipsing X-ray system (Prestwich et al. 
2007), therefore geometric arguments imply that the black hole 
would be eclipsed for i > 78° for the WR properties derived 
by Clark & Crowther (2004). If we adopt a radius of ~ 0.5rL 
for the accretion disk, an eclipse of the X-ray emitting accretion 
disk would require i > 80° . NGC 300 X-1 does exhibit signif- 
icant X-ray variability, but lacks a deep X-ray eclipse (Carpano 
et al. 2007b). Therefore, geometric arguments appear to rule out 
inclinations that would cause a total eclipse of the accretion disk 
(i < 73 ± 2°). However, a glancing eclipse would require i > 
63.5 ± 3.5° for the range of WR radii obtained here. We therefore 
adopt i = 80 - 90° for IC 10 X-1 and i = 60 - 75° for NGC 300 
X-1. 

In Fig.[6]we present the host galaxy metallicity as a function of 
black hole masses for NGC 300 X-1 and IC 10 X-1, plus those for 
all HMXB systems for which the presence of a black hole is unam- 
biguous, whose companion is an OB star with mass M2 > 5Mq, 
i.e. LMC X-3 (Val-Baker et al. 2007), LMC X-1 (Orosz, Steeghs 
& McClintock 2009), M33 X-7 (Orosz et al. 2007), V4641 Sgr 
(Orosz et al. 2001) and Cyg X-1 (Gies et al. 2003). We limit our 
sample to classical HMXB, to ensure that their present-day (oxy- 
gen) metallicities are consistent with the formation of these black 
hole binaries. The majority of black holes in Low Mass X-ray Bi- 
nary (LMXB) systems have masses close to WMq (Remillard & 
McClintock 2006). 

It may be significant that both Wolf-Rayet/black hole systems 
are located in metal-poor galaxies. IC 10 has an oxygen content 
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Figure 6. Comparison between inferred compact object masses. Mi , versus 
metallicity for all HMXB with Mi > 3Mq and M2 > 5Mq. Black hole 
masses infeiTed for NGC 300 X-1 (IC 10 X-1) relate to a WR mass of 2 1 1| 
Mq (25tg^ Mq) and an orbital inchnation of 60-75° (80-90°). 



of log (O/H) + 12 = 8.1 (Gamett 1990) while NGC 300 X-l/#41 
is located at a de-projected distance of 0.43 po, where po = 9.75 
arcmin (Schild et al. 2003). According to Urbaneja et al. (2005), 
the oxygen content at this galactocentric distance in NGC 300 is 
log (O/H) + 12 ~ 8.44, i.e. relatively similar to the LMC for which 
log (O/H) + 12 ~ 8.37 (Russell & Dopita 1990). The only other 
HMXB whose black hole mass is known to greatly exceed 10 A/0 
is M 33 X-7 (Orosz et al. 2007), for which a near identical oxygen 
content of log (O/H) + 12 = 8.42 is inferred at its location in M 33 
from the calibration of Magrini et al. (2007). 

High black hole masses require that the progenitor star was 
very massive and experienced low mass-loss rates (Belczynski et 
al. 2009). Weak stellar winds is a natural consequence of low metal- 
licity (Mokiem et al. 2007). However, orbital periods of IC 10 X-1 
and NGC 300 X-1 are so short that the radius of the black hole 
progenitor star must have been larger than the present separation of 
the components. As such, the progenitor would have experienced 
extreme mass-loss through Roche-lobe overflow. Therefore, recon- 
ciling high black hole masses with close orbital separations is a 
major challenge for binary evolution models. 

In the standard picture, such systems involve a common- 
envelope phase, which would naturally lead to a merger (Podsi- 
adlowski, Rappaport & Han 2003). Alternatively, de Mink et al. 
(2009) propose that the short orbital period results in tidal-locking 
of the stellar rotation, causing a chemically homogeneous evolution 
through rotational mixing (Maeder 1987). In this scenario, binary 
components would remain compact and so circumvent the high 
mass transfer rates of Roche-lobe overflow systems. 

If NGC 300 X-1 and IC 10 X-1 were to survive their second 
supernova explosion, they would form binary black hole systems, 
merging on a timescale of a few Gyr. Binary black hole mergers 
has been considered by Sadowski et al. (2008), who argued that 
their detection rate may be much higher than double neutron star 
systems for current gravitational wave experiments. 

In conclusion, new VLT/FORS2 time-series spectroscopy of 
the WN star #41 in NGC 300 is presented, which confirm that it 
is physically associated with the NGC 300 X-1 system. We find 
that NGC 300 X-1 hosts the most massive stellar-mass black hole 
known, with the exception of the other extragalactic WR/black hole 
systemic lOX-I. 
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